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1. RELATIVISTIC ORDER PARAMETERS 

Relativistic effects in materials containing heavy atoms pervade many areas of 
condensed-matter physics [1-4]. The present chapter is devoted to the theory of 
relativistic effects in a particularly interesting class of materials, namely supercon- 
ductors and other superfluids [5-9]. At its present stage that theory is limited to 
single-particle effects (on the level of the Dirac equation); a relativistic description 
of the particle-particle interaction is not attempted. However, already on the single- 
particle level several unexpected relativistic phenomena in superconductors emerge, 
among them a novel type of spin-orbit coupling, present in superconductors only. 
The origin and the physics of these phenomena are discussed in this contribution. 

The fundamental ingredient in any description of superconductivity is the order 
parameter (OP). The BCS theory of superconductivity [10], in its original formula- 
tion, is not easily generalized to the relativistic domain. Effects of spatial inhomo- 
geneity of the lattice and the OP, of different OP symmetries, and of magnetic fields 
are, however, easily incorporated within the framework of the Bogolubov-de Gennes 
(BdG) reformulation of the BCS theory [11], which also provides a suitable basis for 
the relativistic generalization. In the conventional BdG theory the OP describing 
the Cooper pairs can be written in terms of products of field operators, ^(r^j^r'), 



ifjl(r)i/j^(r'), -0-|-(r) , i/'|(r / ), and ipi(r)ipi(r'), where the ifj a are components of the two- 
component Pauli spinor ^. (All spatial arguments r are actually vectors r, but for 
ease of notation the arrows are suppressed here and below.) Linear combinations of 
these products can be formed to represent the usual singlet and triplet states, which 
are even or odd under interchange of r and r', respectively [6,11-14]. The singlet pair 
is described, for example, by ^ T \r)ia y ^{r') = ^(r^j^r') — ^(r)^^'), which is 
manifestly even under interchange of r and r', and odd under that of spin up and spin 
down, as expected for a singlet state. Since the Pauli matrix ia y is the nonrelativistic 
time-reversal matrix, this way of writing the pair additionally brings out clearly that 
the BCS Cooper pair consists of time-conjugate single-particle states [11,13]. 

A convenient way of describing arbitrary pair states is by expanding the general 
2x2 pair potential matrix A(r, r') in the set of matrices ia y and iaa y , where a 
is the vector of Pauli matrices, and the three components of iaa y are the matrices 
appearing in the Balian-Werthamer parametrization of triplet OP [15], 

^ T (r)A(r, r')^(r') = A (r, r')i> 7 \r)ia y ^{r') + A(r, r')i> 7 \r)iaa y ^(r'). (1) 

Ordinarily one would use the unit matrix I and the vector of Pauli matrices a as 
basis for this type of expansion, since by using this symmetry adapted set of matrices 
one achieves a separation of scalar quantities 1^) from quantities transforming as 
a vector (^cr^), i.e, a classification with respect to irreducible representations (IR) 
of the rotation group. However, Cooper pairs are formed with two creation (or two 
annihilation) operators, and not with one creation and one annihilation operator, 
and a scalar-vector separation is only achieved in terms of the Balian-Werthamer 
matrices [15]. (^> T a^> , for example, does not transform as a vector under rotations, 
while ^ T iaa y ^> does.) 

In a relativistic formulation the spin index a is replaced by the label for the 
components of four-component Dirac spinors. Out of these four components one 
can form sixteen pairs, which can again be represented in terms of a complete set of 
matrices. One might be tempted to employ the sixteen 4x4 Dirac 7 matrices for this 
purpose, since these usually lead to a separation into so called 'bilinear covariants', 
i.e., objects transforming under some IR of the Lorentz group C [16]. However, just 
as in the nonrelativistic case, the usual set of matrices does not achieve this if the 
object to be expanded is formed with two field operators of the same kind. A set of 
matrices in terms of which such a separation into IR of C is achieved, was derived in 
Ref. [6] , where it was shown that the usual classification into five bilinear covariants, a 
Lorentz scalar (one component), a four vector (four components), a pseudo scalar (one 
component), an axial four vector (four components), and an antisymmetric tensor 
of rank two (six independent components), carries over to the superconducting OP 
if the Dirac 7 matrices are replaced by a new set of matrices, denoted f\. Explicit 
expressions for these f) matrices are given in Ref. [6]. These five OP, with a total 
of sixteen components, exhaust the possible pairings which can be formed from two 
Dirac spinors. As an example, consider the most important OP of all, the BCS singlet 
OP, which can be written in terms of Pauli spinors as ^ BC s,non-rei _ ^ T j(j y ty its 

relativistic generalization reads ^ BCS > rel = ^) T f)Q^^ where t)q = ■y 1 ^ 3 is one of the 
sixteen fj matrices (the one yielding a Lorentz scalar) , the \& are now four-component 
Dirac operators, and 7 1 and 7 s are two of the standard Dirac matrices. 



Initially [5], the relativistic generalization of the BCS singlet OP was constructed 
by replacing the nonrelativistic time-reversal matrix t = ia y in ^ BC s,non-rei 
relativistic one, T = 7 1 7 3 , when constructing ^ BCS > rel 5 i. e . 7 it was postulated that 
relativistic Cooper pairs still consist of time-conjugate states. In later work [6], this 
turned out to be unnecessary, since the matrix rjo = T automatically emerges as that 
of the sixteen fj matrices leading to covariant pairs, which reduces to the BCS singlet 
OP in the nonrelativistic limit. It is important to point out that the relativistic 
generalization of the BCS OP is uniquely determined by both prescriptions, and that 
both indeed lead to the same result, ^ BCS > rel = ^ T ?y #. 



2. DIRAC EQUATION FOR SUPERCONDUCTORS 



In terms of the relativistic OP one can now generalize the entire theory of su- 
perconductivity to the relativistic domain. At the heart of this generalization are the 
Dirac-Bogolubov-de Gennes equations (DBdGE), which for a BCS-like spin-singlet 
superconductor read [5,6] 

7 °[c7 • p + mc 2 (l - 7°) + q^A^u n (r) + / d 3 r'A(r, r')r] v n (r') = E n u n {r) (2) 



-7°[c7-p + mc 2 (l -7°) + q^A fl ]*v n (r) - J d 3 r' A* {r,r')fjou n {r') = E n v n (r), (3) 

where both u n and v n are four-component (Dirac) spinors, representing particle and 
hole amplitudes, while 7^ is the four vector of 7 matrices (in standard notation, 
and with a summation over repeated greek indices implied), and 7 the corresponding 
three vector, containing 7 1 ,7 2 , and 7 3 . is the four potential, and A the pair 
potential. Equations (2) and (3) generalize the conventional Bogolubov-de Gennes 
equations (BdGE) for spin-singlet superconductors [11] to the relativistic domain [5,6] 
in the same way in which the conventional Dirac equation generalizes Schrodinger's 
equation [3,16]. The potentials A(r, r') and A^{r) in Eqs. (2) and (3) are to be 
regarded as effective potentials incorporating the electron-electron interaction in the 
spirit of mean-field [11] or density- functional [12] theory. 

The pairing of time-conjugate pairs, described by the matrix t\q, is, in fact, char- 
acteristic only of BCS-like spin-singlet superconductors. For triplet superconductors, 
even in the nonrelativistic case, one needs parity conjugation p in addition to time 
reversal t, in order to express the order parameter and the underlying pairing state in 
terms of fundamental discrete symmetries. From the two operations, p and t, acting 
on two-component (Pauli) spinors one can form four antisymmetric pairing states, 
representing singlets and triplets [6,13]. The relativistic generalization, P and T, of 
the symmetry operations p and t is, of course, well known [3,16] and the above non- 
relativistic prescription for forming pair states is easily generalized to the relativistic 
case [6]. Apart from P and T, relativity also provides a third fundamental discrete 
symmetry, charge conjugation C, which matter of principle, also be used to 

form pair states. From the three operations P, T, and C, acting on four-component 
Dirac spinors, one can form sixteen antisymmetric pair states, which are nothing 



p 



Figure 1. Solid lines: Schematic energy spectrum for a relativistic BCS-like 
superconductor with a point-contact pair potential, in a constant external potential 
(i.e., disregarding lattice and band structure effects), calculated by analytical diag- 
onalization of the DBdGE (2) and (3). In order to make the superconducting gap 
and the relativistic gap visible on the same scale, A and fi were artificially enlarged. 
Dashed lines: Conventional nonrelativistic BCS energy spectrum for the same values 
of A and fi. More details are given in the main text. 

but the sixteen combinations entering the five superconducting bilinear covariants, 
obtained independently above, from considerations of Lorentz covariance. 

The Dirac-Bogolubov-de Gennes equations incorporating all these OP read [17] 

15 „ 

f{c^-p + mc 2 {l-f) + qYA^u n {r)-Y^ / d 3 r'A i (r 1 r')ffjv n (r') = E n u n (r) (4) 

i=0 J 
15 „ 

-f[c^f-p+mc 2 (l-f)+q^A fl ]*v n (r)+J2 / d 3 r'A^(r,r')fjfv n (r , ) = E n v n (r), (5) 

i=o J 

where the sums are over all sixteen matrices fji , and the corresponding pair potentials 
Aj. These equations thus contain all five possible OP compatible with Lorentz co- 
variance. In practice one expects any given interaction leading to superconductivity 
to produce only one of these covariant OP. The DBdGE for each of these are easily 
obtained from Eqs. (4) and (5) by setting all pair potentials A;, except for those 
belonging to the IR under consideration, equal to zero. (In the case of the Lorentz 
scalar one recovers Eqs. (2) and (3) in this way.) 

Eqs. (2) and (3), can be diagonalized analytically [5] for spatially uniform su- 
perconductors, employing a point-contact pair potential A(r,r') = Aod(r — r') and 



setting A M = (v, 0, 0, 0), where Aq and v are constant in space. The resulting energy 
spectrum is displayed in Fig. 1, together with its nonrelativistic limit, the conven- 
tional BCS spectrum for the same choice of A(r, r') and \i. For small momenta and 
energies the relativistic and nonrelativistic spectra are very similar, but the relativis- 
ts superconducting gap is slightly shifted to higher momenta with respect to the 
nonrelativistic one. An analytical expression for this shift was obtained in Ref. [5]. 
For larger momenta and energies the relativistic spectrum aquires features of a free 
particle solution of the Dirac equation, such as a linear dispersion relation. 



3. PAULI EQUATION FOR SUPERCONDUCTORS 

Although the task to set up a relativistic generalization of the Bogolubov- 
de Gennes formulation of the BCS theory is formally completed with the derivation 
of the Dirac-Bogolubov-de Gennes equations and all possible order parameters they 
may contain, the underlying physics comes out much clearer by leaving the level of 
the Dirac equation and proceeding to the weakly relativistic and the nonrelativistic 
limits. These limits are important because (i) the nonrelativistic limit of the DBdGE 
must be the conventional BdGE, hence recovering this limit constitutes an essential 
test of the theory, and (ii) in the weakly relativistic limit (defined as including terms 
up to second order in v/c, where v is a typical particle velocity) one obtains the first 
relativistic corrections to the conventional theory of superconductivity. 

By systematic elimination of the lower components of the DBdGE (i.e., those 
which are suppressed by factors of v/c in the weakly relativistic limit) one obtains 
from Eqs. (2) and (3) the pair of equations [5,7] 

[h(r) + 5h(r)} u n {r) + J d 3 r' [A(r, r')ia y + 8A(r, r')} v n (r') = E n u n (r) (6) 

- [h(r) + 5h(r)}* v n (r) - J d 3 r' [A*(r',r)ia y - 6A^r',r)] u n (r>) = E n v n (r) (7) 

where u n and v n are two-component (Pauli) spinors, h(r) = -^[p — f^.( r )] 2 + v(r) — 
fiaB(r) is the normal-state Hamiltonian, including the vector potential A and the 
magnetic field B, and fi = Tiqj (2mc) is the Bohr magneton. Neglecting the relativistic 
correction terms 5h and 5 A, one obtains from Eqs. (6) and (7) the traditional BdGE 
for spin-singlet superconductors [11]. The term 



Sh(r) = 



1 



4m 2 c 2 



h v 4 ' 

ha ■ [Vu(r)] x p-\ V 2 v(r) - — 

2 



(8) 



represents the second-order relativistic corrections appearing already in the normal 
state, i.e., spin-orbit coupling, which involves the gradient of the lattice potential, 
the Darwin term, containing second derivatives of the lattice potential, and the mass- 
velocity correction, which is quartic in p. Similarly, 



5A(r,r') = 



1 



4m 2 c 2 



ti 2 9 
ha ■ [VA(r, r')} x $ + — (V + V) A(r, r') 



(9) 



contains the second-order relativistic corrections involving the pair potential, i.e., 
those which appear only in the superconducting state. (The prime on p and V denotes 
a derivative with repect to the primed coordinate.) These terms depend on the pair 
potential in a similar way in which those of Sh depend on the lattice potential, and 
are referred to as the anomalous spin-orbit coupling (ASOC) and anomalous Darwin 
terms, respectively. 

More correction terms arise if additionally triplet OP are included. In this case 
already in the nonrelativistic (zero order in 1/c) terms in Eqs. (5) and (6) one must 
add to the term A(r, r')ia y the sum Ylj= x y z A?'( r ' r ')®j®y [14,15]. Furthermore, 
Eqs. (8) and (9) stem from Eqs. (2) and (3), and hence do not include OP involving 
charge conjugation, neither for singlets nor for triplets. The generalization of Eqs. (8) 
and (9) involving all possible OP discussed in the previous sections is found from 
Eqs. (4) and (5) [17] and reads 



8h(r) = 



4m 2 c 2 



\ vMr) + KS.(V v )*P-l- + ^VuV il 



(10) 



and 



<?A(r,r') = -L [<? (j?- lA{rj) V 2l +V 12 r (f - -/{r'))*' 



+ 



4m 2 c 2 



apD 22 a*p* + — ( V z Pn + DnV 



7 /2 



(11) 



where the 4x4 matrix V is defined in terms of the complete set of fj matrices as 

15 



(12) 



i=0 



and the action of the product T>yiD 21 on a function f(r) should be interpreted as 

15 . 

A2V* 21 f(r)= J2 / dr ,/ A i (r // ,r)AJ(r / ,r // )/(r , )»/} 12 »7Si- (13) 

i,j=0 J 



The subscripts in, e.g., T>i 2 , refer to 2 x 2 sub matrices, i.e., T>\ 2 is the upper right 
block of T>. If all pair potentials but that multiplying 170 are se t equal to zero (i.e., 
if one specializes to a BCS-like spin-singlet OP) these equations reduce to (8) and 
(9). On the other hand, through the OP matrices f}\ . . .1715 and the corresponding 
pair potentials they contain all relativistic corrections for singlet and triplet super- 
conductors up to second order in 1/c, including those forming their Cooper pairs via 
charge conjugation. 

We now turn to a discussion of the physics behind the weakly relativistic cor- 
rection terms of Eqs. (8), (9), (10), and (11). We start by considering only the cor- 
rections appearing for BCS-like singlet superconductors, i.e., those given in Eqs. (8) 
and (9). The spin-orbit, mass- velocity, and Darwin corrections to h(r), contained 
in Eq. (8), appear already in the normal state and are well known from other ar- 
eas of physics [3,4,16]. Although these corrections are thus not of superconducting 
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Figure 2. Difference in power absorption AP S between left-hand circularly 
polarized light and right-hand circularly polarized light in the superconducting state, 
divided by the corresponding normal-state difference AP N , as a function of temper- 
ature divided by the transition temperature T c . The difference AP S is proportional 
to $t[frjj! (uj, T, H)], the imaginary part of the offdiagonal elements of the conductivity 
tensor in the superconductor, as a function of frequency, temperature and magnetic 
field (oj = 4.5meV and H = 0.05T for the data in the figure) [19]. 



origin, their effect on observables in a superconductor can be dramatically modified 
by superconducting coherence. The conventional spin-orbit coupling (SOC) term in 
Eq. (8), for example, breaks rotational symmetry of the spin degrees of freedom with 
respect to the orbital degrees of freedom. It is well known that this broken sym- 
metry gives rise to circular dichroism in the magneto-optical response of metals [4]. 
Below the critical temperature T c the dichroic response of a superconductor is dras- 
tically modified. A microscopic theory of this modification was recently worked out 
in Refs. [18] and [19] on the basis of Eqs. (6) to (9), and the results can qualitatively 
account for experimental observations. As a representative result we present in Fig. 2 
a graph displaying the circular dichroic response of a simple model superconductor as 
a function of temperature, normalized by dividing by the corresponding normal-state 
result [19]. Without SOC both numerator and denominator would be zero, while 
without superconducting coherence the curve would be flat throughout. The strong 
peak seen immediately below the transition temperature thus arises only from the 
simultaneous presence of superconducting coherence and SOC [19]. 

The mass-velocity term, on the other hand, provides a small correction to the 
bare electron mass. It attests to the ingenuity of experimental physicists that the 
very small change of the Cooper pair mass in a superconductor due to the relativistic 
mass enhancement has already been measured [20,21,22]. Another consequence of 
this term in superconductors is the small shift of the energy gap of a homogeneous 



superconductor, mentioned at the end of Section 2. 

The anomalous spin-orbit and Darwin corrections to A(r, r'), given in Eq. (9), 
are fundamentally different from those to 8h(r), since they depend explicitly on the 
pair potential and are nonzero only in the superconducting state of matter. These 
terms were derived for the first time in Ref. [5]. The anomalous spin-orbit coupling 
provides a contribution to dichroism in superconductors, which can be distinguished 
from that of conventional SOC due to their very different temperature dependence 
[18,19]. The first appearence of an SOC term containing the pair potential, i.e., the 
term nowadays called ASOC, dates back to 1985, when Ueda and Rice postulated such 
a term on group theoretical grounds in their phenomenological treatment of p-wave 
superconductivity [23]. However, at that time it was not clear how such a term could 
be obtained microscopically, and what its detailed form was. These questions were 
answered only ten years later on the basis of the theory outlined above [5]. Both 
Darwin terms, the conventional and the anomalous one, have also been rederived 
phenomenologically [9]. This rederivation showed that the anomalous Darwin term 
can be understood as a consequence of relativistic fluctuations of paired particles in 
the pair potential of the superconductor, in a similar way in which the conventional 
Darwin term can be understood as a consequence of fluctuations of charged particles 
in the electric (lattice) potential [16]. 

We now turn to the relativistic corrections appearing for more exotic supercon- 
ductors, contained in Eqs. (10) and (11). A closer look at these equations reveals that 
there are several terms which are straightforward generalizations of those in Eqs. (8) 
and (9), but also two terms which do not have counterparts in those equations. The 
first of these is the term containing the pair potential operator in Sh(r). It is a new 
feature of these equations, not present in other Bogolubov-de Gennes equations, that 
terms related to superconductivity appear also on the diagonal of the BdG matrix. 
The second is the cross term in 5A(r, r'), containing both the pair potential operator 
and the vector potential. Although this latter term has the usual minimal coupling 
form (in which p is replaced by p — (q/c)A) it is, in fact, highly unusual, since it 
leads to a term containing the pair potential multiplied by the vector potential. It 
is of considerable fundamental interest that such unusual coupling becomes possible 
in the theory of superconductivity, but detailed consequences of this remain to be 
explored. Both of these novel terms contain the offdiagonal blocks of the matrix 
"D, which are nonzero only for OP formed with relativistic charge conjugation as 
fundamental discrete symmetry. These OP do vanish in the strictly nonrelativistic 
limit [6,7,17], and are thus intrinsically relativistic. Interestingly, for these OP there 
are no anomalous Darwin and spin-orbit terms, since these involve only the diagonal 
blocks of the matrix X?, which are zero for the OP involving charge conjugation. In 
the absence of this purely relativistic type of pairing the matrix T> is block diago- 
nal, and all terms containing T>i2 and T>2i vanish. The remaining terms provide the 
leading order relativistic corrections to the nonrelativistic singlet [11] and triplet [14] 
Bogolubov-de Gennes equations. 

To conclude this section we emphasize that all results up to this point are in- 
dependent of the form of the interaction leading to superconductivity. The role of 
the interaction is to select among the possible OP, discussed in Section 1, the one 
which is realized in a given superconductor. Relativity alone does not determine what 
interactions can give rise to superconductivity, but it puts strong constraints on the 



form and transformation behaviour of the resulting OP (e.g., it limits their number 
to sixteen, forming the five covariant combinations listed above). Furthermore, re- 
gardless of the detailed nature of the interaction, if there is superconductivity at all, 
then there will be relativistic corrections to it, and the form of these is universal, i.e, 
independent of the interaction itself. For example, any superconductor with a BCS- 
like singlet OP, no matter what (quasi)particles are paired and by what interaction, 
must display the relativistic corrections contained in Eqs. (8) and (9). 

4. EXPERIMENTAL ASPECTS AND OUTLOOK 

By now there are several predictions for experiment, which have been extracted 
from the general theory outlined above. These are (i) the existence of new types 
of order parameters in a relativistic theory, which vanish in the nonrelativistic limit 
[6,7,17], (ii) a small shift in the position of the energy gap [5], (iii) relativistic fluc- 
tuations in the pair potential, giving rise to the anomalous Darwin term [9], (iv) a 
drastic modification of the contribution of SOC to the dichroic response of supercon- 
ductors below T c [18,19], (v) a contribution of ASOC to the same dichroic response 
[18,19], (vi) the influence of relativity on the OP symmetry, even for conventional 
superconductors [5,6,7], and (vii) the relativistic change of the Cooper pair mass 
[5,7,20,21,22]. Clearly, the task to develop the theory to the point at which mean- 
ingful comparison with experiment can be made is not an easy one, and only first 
(preliminary, but encouraging) steps could be reported in this chapter. However, a 
few general considerations may facilitate experimental detection of these effects: 

First of all, since relativistic effects become important at high velocities, one 
should study superconductors in which the electrons move rapidly. This is the case 
either if the normal-state metal has a very high Fermi velocity, or if there are suf- 
ficiently heavy elements in the lattice, whose orbitals are hybridized with the con- 
duction band. Examples of such systems are the heavy- fermion superconductors like 
UPt^ - containing some of the heaviest elements found in solid-state compounds, the 
high-temperature superconductors - which include between the CuO layers heavy el- 
ements like Ba, La or Hg, or BCS superconductors with heavy atoms in the lattice, 
such as Pb. Apart from high velocities, rapid spatial variations of the pair poten- 
tial particularly favour the anomalous relativistic effects (ASOC and the anomalous 
Darwin term), since these depend on derivatives of the pair potential. Situations 
in which the pair potential varies rapidly are, e.g., superconductivity in thin films 
and small grains, multilayers, and the vortex state. In all these cases a short coher- 
ence length allows more rapid spatial variations, and thus increases the anomalous 
relativistic effects (an observation which points again at high-temperature supercon- 
ductors). Presumably, the most important relativistic corrections in these materials 
are those of Eq. (8) (which are already routinely taken into account in relativistic 
band-structure codes for the normal state, but not included in most approaches to 
superconductivity). However, new and interesting phenomena can be expected to 
arise as a consequence of the anomalous corrections of Eq. (9). 

A more detailed exploration of the fascinating territory of the interplay of su- 
perconducting coherence and relativistic covariance remains a task for the future, 
but already now it seems certain that the two areas can benefit from each other 



both in the study of down-to-earth superconductors, which was the subject of this 
paper, and in the field of cosmological and high-energy analogies of coherence and 
superconductivity [24], to which many of our results also apply. 
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